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K — >■ 7T7T Electroweak Penguins in the Chiral Limit 
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We report on dispersive and finite energy sum rule analyses of the electroweak penguin matrix elements 
((7T7r)2|Q7,8|^°) in the chiral limit. We accomplish the correct perturbative matching (scale and scheme de- 
pendence) at NLO in a s , and we describe two different strategies for numerical evaluation. 
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1. Introduction 

The subject of this talk is the calculation of 
the K — > 7T7T matrix elements of the electroweak 
penguin operators, conventionally denoted as Qr 
and Qg Jy]. Interest in these quantities is twofold. 
First, the Standard Model prediction of e'/e de- 
pends crucially on the K-to-27r matrix elements 
of the QCD penguin (Q 6 ) and EW penguin (Q 8 ) 
operators, @,|,|: 
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where the above matrix elements are evaluated 
in the MS-NDR renormalization scheme at scale 
fi = 2 GeV. Secondly, the leading term in the 
chiral expansion of (Q7,s) can be related to the 
V — A QCD correlator with flavor ud, whose spec- 
tral function is experimentally accessible through 
r decays This provides the opportunity 

to perform a data-driven evaluation of hadronic 
matrix elements, whose results can then be con- 
fronted with other non-perturbative techniques 
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(like Lattice QCD and 1/N C expansion). There- 
fore electroweak penguins provide an interesting 
theoretical laboratory for kaon physics. 

2. Summary of formal results 

Our first goal is to relate the K — > tttt matrix el- 
ements to the observable spectral functions. The 
first large step is done for us by chiral symmetry, 
which relates {t^^Qt ,s\K) to vacuum expectation 
values of two 4-quark operators Q 

2 (Oi) , 



{(mr) 2 \Q 7 \K°) p o = 
<(7T7r) 2 |Q 8 |^ > p o = 
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where Fq is the pion decay constant in the chi- 
ral limit and the definitions of 0\ and 0$ can 
be found in Ref. ||. The above relations receive 
NLO chiral corrections (whose order of magni- 
tude is ~ m^-/(4rrF w ) 2 = 0.18) which have been 
investigated within chiral perturbation theory . 

The second step consists in relating the v.e.v.'s 
of Oi, 8 to the V - A QCD correlator AU ud (s) = 



Ily 1 \s) — n^ l; (s), with flavor structure udH 
The outcome is that: 

• (Oi) is related to the integral of s 2 AII(s) 
over all spacelike momenta; 
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• (0$) is related to the leading singularity 
of An(s) at short distance (large spacelike 
momenta) . 

Essential tools in a formal derivation of the above 
results and in relating (Oi,s) to the spectral func- 
tion are the dispersive representation of the cor- 
relator 



Q 2 



and the operator product expansion (OPE) rep- 
resentation, valid for large space-like momenta 
Q 2 3> Aq CD , which, through 0(a 2 ), has the form 



An 



OPE 



(Q 2 ) 
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d>2 



Q2 

ad(p) + bd{n) In -^7 
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a ( i{p.) and bd(p) are combinations of vacuum ex- 
pectation values of local operators of dimension 
d. In the chiral limit, the above sum begins with 
0,6, whose NLO expression is 



a 6 (p) = {2na s + A 8 a 2 s ){O s ) + A ia 2 s {Oi) 
The formal analysis is summarized by 



M(a s 



h(p)-H s (p) 
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where the matrix M(a s ) contains the matching 
coefficients at NLO in a s , and encodes the correct 
renormalization scheme dependence. Its explicit 
form can be found in Refs. (8j. 7x,8 are related to 
the spectral function as follows, 
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while H l s are related to subleading (higher di- 
mensional) terms in the short distance expansion: 



Hi(p) 
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5 The coefficients Ai t g depend on the definition used for 
75 in d ^ 4. They can be found in Ref. fa]. 



A close look at the input needed in Eq. (||) and 
at the database on Ap (Refs. j|,D) shows that 
data alone do not provide enough information for 
a reliable determination of the spectral integrals. 
In fact, we run out of data at s = m 2 , while the 
weighs entering 7x,8 are growing with s. We there- 
fore need to supply some extra theoretical input 
to optimize the use of data in this problem. Along 
this path, we have followed two approaches, which 
can be schematically characterized as follows: 

I. Use of data plus QCD integral constraints 
(classical chiral sum rules). |8j 

II. Use of data plus QCD duality (Finite En- 
ergy Sum Rules, FESR). § 

3. Numerics I: Residual Weight Method 

This approach is an efficient way to enforce the 
constraints from classical chiral sum rules 
in the evaluation of the spectral integrals Ii^(p), 
which have the generic form: 



I{p) = / dsK(s,p)Ap(s) . 



(9) 



The basic idea is to write K(s,p) as the sum of 
an arbitrary combination of the weights 1, s, and 
slogs and a residual term r, as follows: 

K(s,p) = xl + ys + zs log s + r {x ^ y ^ z} (s,p) .(10) 

From Eq. ( pij| ) and the chiral constraints one then 
has 



I(p) = F 2 (x-z 
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dsr{s,p) Ap(s),(ll) 



and can choose (x,y,z) (so far arbitrary) in such 
a way as to minimize the total uncertainty (asso- 
ciated with both the theoretical input for F and 
5M 2 , and the large errors on Ap(s) for s > 2.5 
GeV 2 ). The method is effective in those situa- 
tions where choices for x,y, z exist which make 
r ~ in the region with no data (but before the 
onset of asymptotia) without amplifying the ef- 
fect of imprecise knowledge of Fq and SM 2 . Pro- 
cedural details and intermediate results can be 
found in Ref. Here let us recall that the 

method leads to a reasonably accurate determi- 
nation Jx 8 at different values of p between 2 and 
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4 GeV (with accuracy deteriorating at large fi). 
This method tells us very little about the sub- 
leading terms -ffi,8, which can be neglected only 
for sufficiently large /x. We have assumed that 
#i, 8 (>t = 4GeV) = in Eq. (|) and used the 
NLO renormalization group equations to evolve 
the results down to the lower scale fx = 2 GeV. 
Numerical results are reported in Sect. || 

4. Numerics II: FESR 

In principle, this method allows us to determine 
the coefficients ad appearing in the OPE expan- 
sion for the correlator AII(s) in the asymptotic 
Euclidean region (Eq. (||)). The determination of 
0,6, in particular, provides a direct extraction of 
(0$) at [i = 2 GeV. The determination of the ad, 
d > 8, also allows one to estimate the combina- 
tions #1,80" = 2 GeV), needed as input to the 
lower scale version of the RWM analysis, Eq. (Q). 
Let us now summarize our implementation of the 
FESR analysis @. 

An FESR analysis, applied to a generic corre- 
lator II(s) with spectral function p(s) is based on 
the following relation (consequence of Cauchy's 
theorem, and the analytic structure of the corre- 
lator) : 



J[w, so] + R[w, s ] = fw[so> Wd}} 



Eq. (y_2|) involves the spectral integral 

/•so 

J[w, s ] = / ds w(s)p(s) , 



(12) 



(13) 



the integrated OPE (function of the various con- 
densates) 



fw [s , {a d }] = - 



2ni 



s\=so 



ds w(s)n OPE (s) ,(14) 



and the remainder term (arising from the differ- 
ence n(s) — IIop£;(s) on the \s\ = s circle) 



R[w, s ] 



1 

27ri 



ds w(s) (n c 
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Eq. (|l2] ) is valid for any value of so and any poly- 
nomial weight. However, the presence of R[sq, w] 
pollutes the extraction of the OPE coefficients 
ad in terms of data through Eq. ([l^). There- 
fore it is highly desirable to work with a range 



of so values and with weights w(s) such that 
R[w,sq] « J[w,sq]. To accomplish this, we 
rely on the observation that for |s| large enough 
(s 3> Aq C£ >), the OPE provides a good represen- 
tation of the full correlator along the whole circle, 
except in a region localized around the time-like 
axis. The physics of this breakdown is given by 
the arguments of Poggio, Quinn and Weinberg 
p2j] . As a consequence one expects that weights 
with a zero at s — sq, de-emphasizing the region 
where the OPE fails, are good candidates to gen- 
erate a small-sized R[w, sq]/J[w, sq] B. We use 
this as the basic guiding principle in the selection 
of weights for our analysis. The weights were also 
chosen to avoid large cancellations in the spectral 
integrals, which would lead to a poor statistical 
signal. 

The lowest degree weights satisfying our cri- 
teria have degree 3. Two useful cases (with 
V = s/s ) are: 



My) = (i - vf (i - 
Mv) = i l -y) 2 y 



(16) 
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These choices imply that f Wl2 [sq, {ad}} depend 
on the dimension six and eight condensates, 
thereby allowing us to extract information on a§ 
and as from a least-square fit. The window of so 
values used in the fit procedure is 1.95 GeV 2 < 
so < 3.15 GeV 2 . The lower endpoint has been 
determined by trying ever lower values of so until 
the extracted ae ceases to be consistent with the 
values obtained by the smaller analysis windows. 
From the ALEPH data we obtain (at fi = 2 GeV) : 



a 6 = (-4.45 ± 0.7) • 1CT 3 GeV 6 
a s = (-6.2 ± 3.2) • 1CT 3 GeV 8 



(18) 
(19) 



with a correlation coefficient c(a6,a$) — —0.99. 
Although no precise answer emerges for as, we 
see that a§ (directly related to Qg) can be de- 
termined fairly well. The errors reported above 
are essentially of statistical nature. In order 
to address potential systematic effects due to 
R[sq,w] ^0 (duality violation), we have repeated 
the fit procedure using non-overlapping so sub- 
windows. We find that the fitted parameters in 
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rting evidence for this fact has been found in 
J, within an analysis of the V and A correlators. 
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the different analyses are consistent, thereby con- 
firming that the effect of R[sq, w] is suppressed in 
this case. A more explicit portrait of the FESR 
machinery is obtained by plotting J[wi t 2, so] and 
fwi 2 i s o', £26; a &] as a function of So, as in Figs. |^J|. 
The excellent match of the OPE curve versus data 
increases our confidence that in this analysis du- 
ality violation is under control. 




Figure 1. f Wl [sQ,ae,ag] (continuous curve) and 
J[wi,so] (data points) versus sq (GeV 2 ). 
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Fi gure 2. f W2 [so, osg, og] (continuous curve) and 
J[ w 2> s o] (data points) versus s (GeV 2 ). 



Results for ((7T7r) /=2 |Q 7 , 8 |X ) 



Having outlined the two numerical strategies 
in the above sections, in Table [l] we report (0,7,%) 
in the MS-NDR renormalization scheme at scale 
H = 2 GeV, from methods I, II, and VSA 0. 

Procedures I and II lead to consistent results, 
with II having reduced uncertainty. The underly- 
ing reason is that method II exploits best the low- 
error part of experimental data. Both our results 
seem to point to moderate violation of VSA. For 
a comparison with other recent determinations 
@Mra0Jll see Refs. pj§. Implications of 



Table 1 

MS-NDR results at \i = 2 GeV and n f = 4. 



Method 


(Q 7 )/= 2 /GeV 3 


(Q 8 )/=2/GeV a 


I 


0.16 ± 0.10 


2.22 ± 0.70 


II (ALEPH) 


0.23 ± 0.05 


1.41 ± 0.27 


II (OPAL) 


0.21 ± 0.05 


1.72 ± 0.32 


VSA 


0.32 


0.94 



our results for the value of e'/e in the Standard 
Model will be discussed in B. 
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